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CONVEXITY PROPERTIES OF DISCRETE SCHRODINGER EVOLUTIONS 
AND HARDY’S UNCERTAINTY PRINCIPLE 

AINGERU FERNANDEZ-BERTOLIN 


Abstract. In this paper we give log-convexity properties for solutions to discrete Schrodinger 
equations with different discrete versions of Gaussian decay at two different times. For free 
evolutions, we use complex analysis arguments to derive these properties, while in a perturba¬ 
tive setting we use a preliminar log-convexity result in order to get these properties. Then, by 
proving a Carleman inequality we conclude, in one of the cases under study, a discrete version 
of Hardy’s Uncertainty Principle. 


1. Introduction 


The aim of this paper is to show log-convexity properties for solutions to the free Schrodinger 
equation 

d 

(1) dtUj = iAdUj = i 'y -|- Uj-ep.), 

k=l 

and also for solutions to a perturbed discrete Schrodinger equation 

(2) OppUj — li^AdUj “t" 

when U is a time-dependent bounded potential. 

Our motivation is the relation between log-convexity properties and Hardy’s uncertainty prin¬ 
ciple (see [711^1. 


( 3 ) 


\fix)\ < 1/(01 < with a/3 > i ^ / = 0, 


and in the case a/3 = | then f{x) = Ce “1^1^. This uncertainty principle can be understood as 
an amplification of Heisenberg’s uncertainty principle 

1/2 / /• \ 1/2 


\ 1/2 / ^ \ U2 p 

\xf{x)\‘^dx] / \Vf{x)\‘^dx] > / |/(x)pdx, 

/ XJR'i J JR'* 


where the equality is attained if and only if f{x) = Ce for a > 0. Moreover, writing a 

solution to the Schrodinger equation dtu = iAu as 


i|a:|^/ 4 t , ^ A 


© 


{ityc 

Hardy’s uncertainty principle can be written in terms of solutions to the Schrodinger equation 
in a setting (see M) as follows: 


u(0)||i2(Rd) -F u(l)||i 2 (Rd) < -too, a/3 > — ^ u = 0, 
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SO this result states that a solution to the Schrodinger equation cannot decay too fast at two 
different times simultaneously. The classical proof of this uncertainty principle is based on 
complex analysis arguments (Phragmen-Lindelof’s principle and Liouville’s theorem), but there 
is a series of papers, [HHII] and [5], where the authors prove Hardy’s uncertainty principle in 
this dynamical setting, considering solutions to perturbed Schrodinger equations and using real 
variable arguments. One of the main tools they use is precisely a log-convexity result that 
states that a solution to those equations with Gaussian decay at two different times preserves 
the Gaussian decay at any time in between. This process to prove Hardy’s principle using 
real calculus starts in with a non-sharp result combining the log-convexity property with a 
Carleman inequality, and then in m they use an iterative process to go from this preliminar 
result to the sharp result. 

In the discrete setting, we studied in a previous paper, m, a version of Heisenberg’s principle 
(see also [mi for more references to this uncertainty principle) generated by the discretization 
of the position and momentum operators 

(4) = jhu, = {nh ,... ..., 

k 

where j G lA, = (0, ..., 0 , 1 ,0,..., 0), for fc = 1,..., d, and related it to the discrete 
Schrodinger equation m via a Virial identity. In this case, the minimizer uj (the analogous of 
the Gaussian function) is given in terms of modified Bessel functions of the first kind 

Im{x) = “ y cos(to6») de, m G Z, oj = = {ch\Ij 

Moreover, in the paper we saw how we can recover the Gaussian from the minimizer to 

as the mesh step tends to zero. 

On the other hand, using complex analysis arguments, we gave in m a discrete version of 
Hardy’s principle in one dimension, similar to (H, that can be written in terms of solutions to 
the discrete free Schrodinger equation as in the classical case, 

(5) \uj{0)\ < Ij{a), \uj{l)\ < Ij{P), a-|-/3<2=>u = 0. 



It is natural to think then that we should be able to prove this result following the approach 
in d, so the first thing we need to understand is the analogous version of the log-convexity 
property stated in [5]. Notice that in the discrete setting we can give many discrete versions of 
Gaussian decay. Looking at (H, the first interpretation to play the role of one can think of 

is the inverse of the minimizer ui. However, we can understand the function as the solution 

to the adjoint equation of V/ -I- 2Xxf = 0, the equation satisfied by the Gaussian. If we do the 
same using the operators (H , it is easy to check that now the weight is given in terms of modified 
Bessel functions of the second kind 





e ^ ‘ cosh(mt) dt. 


On the other hand, we can simplify more the discrete interpretation of Gaussian decay, just 
using the weight function . 
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Here we give two different methods in order to prove that solutions to the discrete Schrodinger 
equation with those discrete versions of Gaussian decay satisfy a log-convexity property. Formally 
we see that the log-convexity holds, but trying to justify these formal calculations is where we use 
different methods. First, by relating the discrete solution to a periodic function via Fourier series 
we can use complex analysis arguments in order to justify everything. However, this method is 
only useful when considering solutions to the free case H]). On the other hand, for solutions to ([5|) 
we can give a preliminar log-convexity property, in the spirit of |16j . using a linear exponential 
weight, which, by a simple fact, allows us to prove the log-convexity properties we want in the 
perturbative setting. The advantage of using the first method is that proving the log-convexity 
directly we can get some a priori estimates. These estimates were crucial in the continuous case, 
although in the discrete setting this does not seem to happen. 

Once we have the log-convexity properties for the different weights, since we want to follow the 
approach in [S], we need a Carleman inequality in order to prove Hardy’s uncertainty principle. 
Using modified Bessel functions, it is not clear how we should take the frequency function, but 
in this paper we prove that there is a Carleman inequality associated to for the discrete 

Schrodinger equation, analogous to the inequality in [5], that can be used to prove a non-sharp 
discrete version of Hardy’s uncertainty principle. Hence, we give a discrete uncertainty principle 
whose proof looks like the proof in the continuous case. 

Preparing this manuscript, we learned about a recent and independent result in this direction 
m- There, the authors also prove a sharp analog of Hardy’s uncertainty principle in the discrete 
setting, in terms of solutions to the Id discrete free Schrodinger equation by using complex 
analysis arguments. To avoid the use of complex analysis, and to add a potential to the equation, 
they adapt the log-convexity approach in [S], getting also a non-sharp result in this case. Their 
result improves our version of Hardy’s principle (Theorem 4.1 below), when we focus on the one 
dimensional case, although the proofs of the results are rather different, since here we follow 
more closely the approach in [U] . 

The paper is organized as follows: In Section 2 we give log-convexity properties using the 
weights discussed above for discrete free Schrodinger evolutions, so, in order to justify the formal 
calculations, we use tools of complex analysis. Then, in Section 3, we add a potential and prove a 
result using a linear exponential weight, concluding from this result the log-convexity properties 
of Section 2, now in a perturbative setting. Finally, in Section 4 we prove the discrete version of 
Hardy’s principle by means of a Carleman inequality. 

In our previous papers, we study the discrete Schrodinger equation with the mesh step h, 
so that when h tends to zero the solution to the discrete equation converges to the solution to 
the continuous equation. Here, since we are not going to study convergence of the results when 
h tends to zero, we fixed, just for simplicity, h = 1, although all the results can be written 
introducing this parameter. 


2. Log-convexity properties of discrete free Schrodinger evolutions 

To begin with, we consider that the solution to ([T]) decays when we multiply it by the inverse 
of the discrete minimizer in m- Then we have the following result: 
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Theorem 2.1. Assume u = is a solution to the d-dimensional free Schrodinger equation 

0 which satisfies 

' ' (jJ ■ ' ' OJ ■ 

where uij = Cd,\ 11^=1 ^jk (1/2^) for some A > 0. Then 



In order to prove that F{t) is log-convex in the interval [0,1], we want to use the following 
lemma (see 0). 

Lemma 2.1. Assume that f(t) satisfies dtf = Sf -f Af, where S is a symmetric operator and 
A is a skew-symmetric operator (both independent of t). If [iS,Al] > 0, then H(t) = (/,/) is 
logaritmically convex. 


Proof of Theorem 2.1. Formally, we set fj{f) = and it is easy to check that dtfj = Sfj-\-Afj 
where 




( 333 k-1 



) 93 -ek ) 

\ ^3k 

^3k-lJ ) 





1 E ^ 

- 

' ~ 1 93 -ek 1 

V ^3k 

^3k-lJ J 


where we define, for n G Z, uin = -fn(l/2A). Hence, since F{t) = (/, /), in order to use the lemma 
we need to show that [5, A] = — AS > 0. 


The commutator of these operators is given by 


{SA-AS)g,=lY^ 


^3k^3k+‘2 


U3 


k=l \ \ ‘^3k^3k-\-‘2) 

d / o -2 1.2-4 


Ik-\-l \ I / ^3k^3k-2 

gj-\-2ek + 




+ E 


2i2 


k=l X^^lk + l^L-l 


- - 


<^3k 




^Ik-l 

^3k^3k-2 ^ 

93- 


93-2ek 


Thus, the expression we want to be positive is {[S,A]g,g), which after some calculations is 


y y — ^ — 




^^3k 


EE 

jgZ<i k=l ' 

^3k-23jjk 
' 2^-1 




,^3k + 1^3k-l 


2 ~ ^Jk^ + 




UJ 


3k + l 


+ 


^jk-l^^lk-^-l 
^3k+23bjk 




3k 


CJ. 


3k-l 


“^333 k-2^3 3 k 




\93\- 


Notice that for each k, the expressions that appear multiplying \gj+ek ~ dj-ekl^ and \gj\‘^ are 
exactly the same, so once we prove that this is positive in one dimension, it is straightforward 
to prove it in the general case, so we restrict ourselves to the one dimensional version of the 
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commutator. The first sum is positive by Amos inequality, [U p. 269] or [5H1 (1-9)]: l]{x) — 
> 0 for a; > 0 and j > —1. Notice that since j is an integer, I-j{x) = Ij{x). 
Hence, it remains to prove that the second sum is positive when j G N U {0}. To simplify, we 
will consider x = and divide the expression of the second sum by 4A^. That implies that if 
we are able to prove the following property for modified Bessel functions: 

_ 9,2 , 2 / _ Ij + lix)I,-lix) I,+2{X)I,{X) 

I‘-+i{x)I]_^[x) \lj-i{x)Ij+i{x) Ij{x) 2/2 (a;) 

(6) ^ ^ 

2Ij(x)Ij+2ix) 2/2_^(a:) 2I^_2(.x)Ij{x) j 

for J S N U 0 and a: > 0, then the convexity will hold. 

Since f{x) = x^^'^Ij (a:) satisfies the equation 

we see two behaviors, when x is large enough and x is small enough (with respect to j). Note 
that in the case of the Bessel function of the first kind Jj (x) there is a cancellation term in the 
equation that gives another behavior. 


These behaviors are given by the asymptotics (first for x small enough, then for x large 
enough) 

i^/^y 


Ij{x) 


Ijix) 


y2Trx V 


/ 4j2-1 (4j2 _ i)(4j2 _ 9) (4^-2 _ _ 9)(4^.2 _ 25) 

I ' rvl/o \ O \Q 


8a: 


2! (8a;) 


3!(8a:) 


Therefore, when x is small enough (with respect to j) 

2j + l (2j + 1)(2j3+4j2-3j-4) 


® - 2(2j + l)-x^ 


> 


> 0 , 


(j - l)j(i +1)0' + 2) - (j - I)(j + I)(j + 2) 
if j > 2. The cases j = 0,1 follow a similar argument but the calculations are slightly different. 


On the other hand, if x is large enough the best way to treat m is writing it in the form of a 
quotient of two expressions involving modified Bessel functions. In the denominator we will have 
a product of modihed Bessel functions (which are positive functions). Moreover, since the degree 
in both numerator and denominator is the same, we can avoid the term in the asymptotic 

expansion. 


If we make all the calculations, and only consider the leading term in the expression, we see 
that it behaves like 


2 + 8f 


> 0 . 


Thus, heuristically we have seen that it makes sense to think that m is positive. 

To give a rigorous proof of this, we will use rational bounds for modified Bessel functions in 
order to reduce the positivity of our expression to the positivity of a quotient of two polynomials. 
We will need to treat separately three different cases. For the sake of readability, we avoid the 
calculations here because of the large numbers involved, since the degree of the polynomials 
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involved is quite high. In order to manage all these computations, we use Mathematica as a 
useful tool to substitute coefficients of the polynomials according to the regions we are studying 
in each case. 


First case: j > 17, 0 < a; < In this region we use the following Turanian estimate, 

which is an immediate consequence of [21 Theorem 1]: 

^ ^ < If{x) - Ij-l{x)Ij+i{x) < -^I^{x). 


U + l)\/a:2 + (j + 1/2) 


2-5 


Using these bounds, after some computations we see that the positivity of m depends on the 
positivity of an expression of the following kind, 

(1 + 2j)p{j, a;) - (1 + j)\/l + 4j + ip + ^x‘^q{j, x), 

where p and q are polynomials of positive coefficients whose degrees, as polynomials in x, are 6 
and 5 respectively. Writing the difference (1 + 2j)^p^(j, a;) — (1 + j)^(l + 4j + 4j^ + 4x^)g^(j, x) as 
a polynomial in x (of degree 12) whose coefficients are polynomials in j, we see that in this region 
® is positive. Indeed, we look at the sign of the coefficient of highest degree in x, and when it 
is negative we use x^ < to reduce the degree of the polynomial. Iterating this argument when 
necessary we see that the polynomial is positive. 

Second case: 0 < j < 17, X > 0. Now we use the rational bounds in HZl Theorem 2 and 
Theorem 3] to treat this case. There, the author gives a method to generate upper and lower 
bounds for the ratio based on the completely monotonicity of the function x~^e~^Ij{x). 

More precisely, he proves that there are rational polynomials such that 

(7) 7/j/c^m(x) < r <C Uj,k,m(^x) ^ 

giving also a method to compute these polynomials L and U. Now, setting fc = 5 and m = 0, 
we get a rational bound for (|2|). We can easily check that the denominator of the quotient is 
positive, and, for j < 17, all the coefficients of the polynomial that appears in the numerator (of 
degree 39 in x) of the quotient are positive, so (O is positive in this case. 

Third case: j > 17, x > Here we use again ©, with k = 5, m = 0. As we have pointed 
out in the second case, the denominator of the rational bound is positive, but now when j is 
large some negative coefficients appear in the numerator. Again, the way to prove the positivity 
is to collect the coefficients in x powers and look at the sign of the coefficient of highest degree 
in X. Now the coefficient is positive if j > 17, and using x^ > we reduce the degree of the 
polynomial and start again this process. 

Hence, the separate study of (O in these three cases gives us the positivity of the commutator 
and by Lemma 2.1 we get the log-convexity of F'(t), formally. If we can justify this formal 
argument, then the theorem holds. 


In order to justify all the calculations, what we want to check is that 


OO 

Fc{t)= Y. 

j=-oo 


lujjtr 

I]{a) 


is well defined for t G (0,1), provided that Fq,( 0) -I- Fh(l) < -too and a > 0. Once we prove 
this, we conclude that the formal calculations are valid and then we have that the log-convexity 
holds. 
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As a first step, we prove a weaker result that says that a solution bounded at two times by the 
modified Bessel function is still bounded at any time between them. Then, from this £°° result 
we get the result we want. 


Proposition 2.1. Assume that u = {uj)j is a solution to CD when d = 1 such that 


|wj(0)| + |mj(1)| < CIj{a), Vj e Z, 

for some a > 0. Then there is C > Q such that, for t G (0,1), \uj{t)\ <C{\ 


1 A 1 3 (a) 

i-V yid' 


Proof of Proposition 2.1. We consider that Uj{t) = f{j,t) for a 27r-periodic function /. The 
evolution of f{x,t) is given by 


Moreover, since |/(j,0)| < CIj{a), we have that f{x,0) is extended to an entire and 27r- 
periodic function 

OO 

f{z,0)= z = x + iy, 

j^-OO 

SO f{x^t) inherits these properties for all time. Furthermore, using that 




g-iy _ gO cosh y 


we conclude that \ f{z, 0)| + \f{z, 1)| < for all z = x + iy. Hence 

r sin x sinh y+a cosh y 

\f{x-iy,t)\ < 1 

2(i—1) sin X sinh y+a cosh y 


Since we want something that behaves better than we are going to use, for ?/ > 0, 

the first line when sinx is positive, that is, when x G [0,7r]. On the other hand the second line 
will be useful when sinx is negative, that is, when x G [— tt, 0]. 

We have to distinguish between j positive and j negative, although the procedure we follow 
is the same. The quantity we have to look at is 

f{j,t)=[ f{x,t)e~"^^dx, 

and we have to see that this quantity is controlled by Ij{a). For j positive, we integrate the 
function over the contour described in Figure [2 observing that the integral over the vertical lines 
vanishes due to the periodicity. Thus we see that, thanks to Cauchy’s theorem 


fU,t)= [ f{x-iy,t)e *^)dx, Vy > 0. 

J —TT 


Now we split up the integral in order to use the bounds for |/(x — iy,t)\. Thus, 



CONVEXITY PROPERTIES AND HARDY’S UNCERTAINTY PRINCIPLE 



-§-6-iy 


ly 


Figure 1. Contour integration 


We can write each integral as a difference between a modified Bessel function of the first kind 
and a modified Struve function, both of order zero, having that 

\fij,t)\ < - t) sinhy) - Lo(2(l - t) sinhy) + /o(2tsmhj/) - Lo(2tsmhy)). 

Following the theory in [211 §10.42], we see that 

2 16 
7r(/o(s) - Lo(s)) = - + i?, where |i?| < —. 


Hence, using this expression we obtain that for y large enough. 


^2(1 — t) sinh 2 / 2tsmhy 

On the other hand, from [m Ch. 10, §7], we have 




gO coshy-jy 

2 sinh y 


( 8 ) 


72 ^( 1 + ^ 

gt\/l+aVF-i arcsinh(j7a) 


< 


5j’ 


so, for j large enough we have 


(a) > 


arcsinhO'/a) 

2^J2TTj{\ + 


Thus, if we set y = arcsinh(j /a) (that tends to infinity when j tends to infinity), we have that 
for j large enough. 


la?-j arcsinh(j7a) p3\/^+a^!-j arcsinh(j7a) 

\f{ht)\ < Ct --= Cf 


2jla 


2jla 


^ ^t,oc - 7^5 

VJ 


since v7 and v7(l + behave in the same way as j grows. If j is negative, we use the 

same argument but instead of integrating the function /(x — iy, t) we integrate /(x + iy, t) over 
a similar contour and then we take y = arcsinh(|j|/a). □ 


Since i is not a summable function, we cannot use this proposition directly to justify these 


calculations. Nevertheless, we have that this implies that 

F V <c V 


3^-00 




3^-00 


Ijia + e) 


< + 00 . 
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That the last sum is finite can be seen using the bounds in [5D1 (4)], and the same can be 
done for the time derivatives of Fa+e- Hence we have that for Fa+t the formal calculations are 
correct, so it is a log-convex function for all e > 0. Notice that since the constant in Proposition 
2.1 blows up at t = 0 and t = 1, first we prove the log-convexity of T'a+e in an interval of the 
form C (0,1) and then, by using the convolution expression for the solution. 


Uj{t) = e ^ Um(0)4_m(2if) = e 


_ p-2i(t-l) 


- 1 )), 


m——oo 


m— — oo 


we can let to tend to 0 and ti tend to 1 to conclude the log-convexity in [0,1] In other words, we 
have that 




E 




/2(a-f e) 


< 


< 



Finally, by Fatou’s lemma. 


1^1 ft)I ^ I- l^ijft)!" 


Ij {a 

J^-CO J 


e-J-0 ^^ 

j^-co J 


ft E 


l»j(0)p 

iHa) 


E 


iHa) 


so the theorem holds. The same method can be used to justify the formal calculations in the 
general d-dimensional case. □ 


Now, as we have pointed out in the introduction, we have other interpretations of Gaussian 
decay, so let us consider the solution to the adjoint equation that solves the modified Bessel 
function Ij{x), 

>^jkZj - {Zj+e^ - Zj-ek) = 0, j G 

It is a simple computation to check that now the weight is given in terms of modified Bessel 
functions of the second kind Kj{x). Using this weight, we have the following result: 

Theorem 2.2. Assume u = (ujftgzd is a solution to Q which satisfies 
(9) z^\uj{0)\^ + Y < +00’ 

where Zj = Cd,\ 11^=1 some A > 0. Then 

Hit) = Y is logarithmically convex. 


As before, we are going to prove the log-convexity formally. In order to justify the calculations, 
we can argue in the same fashion as in Theorem 2.1, now proving the following one dimensional 
result (whose proof is based on the same arguments that we have used to prove Proposition 2.1): 


Proposition 2.2. Assume that a solution to the Id discrete Schrodinger eguation ([T]) satisfies, 
Vj G Z, Kjia)i\uji0)\ + |uj(l)|) < C, for some C > 0 and a > 0. Then, for t G (0,1) we have 


Kj{a)\uj{t)\ < C« 


i-d vliT’ 


if j is large enough. 
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Proof of Theorem 2.2. Again, we set fj = ZjUj and carry out all the process to compute [5,-4] 
in this case. We define = Kn{l/2X) for n G Z, noticing that in the previous theorem we have 
done this for the inverse of Zj. Then in this case we have 


fe=i 


2z 


d 

+ EE 


Jk 




4-1 


^3k 




3 k 


l^i+efe do — e 


^3k + l 


k=l 


2^2 

3k 


2^2 

lfc-1 


222 


2%)b%fc+2 


Z3k+-2Z3k 

222 
^^3k + l 


22 

^jkZjk- 


Z3k-2Z3k 

222 

ifc-i 


As before, we only need to prove that the commutator is positive in one dimension. The first 
sum is positive due to the symmetry K-j(x) = Kj{x) for j € N and the inequality K^(x) < 
Kj-i{x)Kj+i{x), valid for j > 0 and a: > 0. 

On the other hand, the positivity of the coefficients in the second sum is not straightforward. 
For simplicity, we define A.j(x) as the j-th coefficient in the second sum, 

K]^,{x)+K]_,{x) K]{x) K]{x) K]^,ix) 

2K]{x) 2Kl,{x) 2K]^,ix)^ 2K,ix)K,+,ix) 

_ Kj+2{x)Kj{x) _ Kj_2{x)Kj{x) 

2K]_,,{x) 2 K,{x)K,_2{x) 2K]_,{x) ’ 

and we need to prove that Aj {x) > 0 for j G N U 0 and a: > 0 . We start proving separately the 
cases j = 0,1 and then we study the case j > 2. 

To prove that Ao(a) > 0 we see that this is equivalent to prove that ^ which is 

a consequence of the estimates given in [J1 Theorem 2]. 

If j = 1, we consider two cases 0 < x < 1.1 and x > 1.1. In the first case we use different 
estimates for each term in Ai(x) 

K2{x) ^ , 3 a:2(x) ^ 1 a:o(x) ^ X 

A:i(x) 2x’ Ki{x)Kz{x) 1 + 4 1 + + 1' 

The first estimate comes from the differentiation of the function c^Kq^x), which is a completely 
monotonic function according to [181 Theorem 5] and the third estimate was proved in [211 
Theorem 1]. Asa result, we obtain that Ai(x) is positive ifp(x) —Vl + 4:x'Zq(x) is positive, where 
p and q are polynomials of degree 6 and 5 in x respectively. As we have done in the previous 
result for the modified Bessel functions Ij{x), studying the sign of p2(a;) — (1 + Ax‘^)q{x)^, which 
is a polynomial of degree 11, we obtain the positivity in the first region. 

In the second region, we consider the following estimates, 

K 2 {x) 8 + f K 2 {x) 8+f Kjjx) 1 

Ki{x) 8-4 + 4’ K^ix) 8+f + f’ Kfix) l + i + ^a’ 

In this case we use the completely monotonicity of e^Ki{x) to prove the first and second 
estimate. The third is given in [S] Theorem 2]. Again, this gives a rational bound for Ai(x) and 
it is easy to see that both numerator and denominator are positive if x > 1.1. In the case of 
the denominator, it can be written as a product of positive polynomials. For the numerator, we 
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obtain a polynomial of degree 11 with some negative coefficients, but using x > 1.1 they can be 
hidden in the positive coefficients. 


If j > 2 we use the recurrence of modified Bessel functions of the second kind 


K,+i(x)-K,_i(x) = 21k,(x) 


to rewrite Aj (x) as 


2 2 ‘^pKHx) 

r%(x)=2/- > 


Kj_,ix)K]^,{x) 


Kj+i{x)Kj_i{x) 

KVx) 


K]{x) 


Kj+i{x)Kj_i{x) 


+ 


_ K,ix)K,_2{x) ^ Kj+,{x) _ K,ix)K,+2ix) \ 


2K,{x)K,_2{x) 2K]_,{x) 2K,{x)K,+2{x) 2K]^^{x) / ' 


As we have done in the case j = 1, we split up a; > 0 in two regions. 
First, if a; > ^, we use the following estimate, given in [31 Theorem 2]: 
1 - Kl{x) , 1 


-r < 

1 + T Ky — lKyj^i 


< 


1+i- 


—, for V > 1/2. 


After using these estimates, we obtain that the positivity of Aj depends on the positivity of 
a polynomial of degree 7 in the variable x. Studying the sign of the coefficient of highest degree 
(which is going to be positive) and using a; > 3j/2 to reduce the degree of the polynomial, we 
obtain that in this region Aj{x) > 0 after some iterations of this argument. Finally, if 0 < a; < 
we change the upper bound, using [m Corollary 1], 


1 ^ K^.jx) ^ 

1 + ~ Ky—lKyj.1 1 -|- 


1 

1 ’ 
V-^ + y/x^ + (v-^)^ 


for V > 1/2. 


Now the positivity of Aj depends on the positivity of an expression of the kind Pj{x) + 
■^1 — 4/ + 4j^ + Ax'^qj{x), where pj and qj are again polynomials of degree 5 and 4 in x. On the 
one hand, using the same strategy as before, we can see that pj and Qj are positive polynomials 
if / > 2, 0 < a; < j. When x > j we compute the difference p^ — {1 — Aj + Aj^ + Ax^)q'j, which 
is a polynomial of degree 9 in x, and we see that this difference is positive (in order to use the 
same argument as before, we have to distinguish j < x < Aj/3 and 4j/3 < x < 3j/2). Thus, 
there is no change of sign in pj(x) + \/l — 4j + 4/^ + Ax'^qj(x) and this implies that Aj(x) > 0 
as we want. 

This completes the proof of Aj{x) > 0, Vj G N U 0, a; > 0, and therefore we have that 
the commutator of the operators S and A is positive, giving as a result the log-convexity of the 
desired quantity provided that all the quantities involved are finite, by using Proposition 2.2. □ 


We can simplify more the discrete interpretation of the Gaussian decay, and use the weight 
function , having the following result: 

Theorem 2.3. Assume u = (wj/jgzd is a solution to the equation (|T]) which satisfies 

(10) ^ + E < +«)> 

jGZ'^ 
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for some A > 0. Then 

G{t) = \uj{t)\‘^ is logarithmically convex. 


Proof. Formally, we consider fj = Uj and compute the operators S and A so that we can 
apply Lemma 2.1. In this case 
( 11 ) 


{[•S,A]f,f) = sinh(2A) Y l/i+fifc - + 2sinh(2A) Y - O' 

jgzd k=i fe=i 


In order to justify the formal calculations we need again an £°° result in one dimension 
analogous to Proposition 2.1 and Proposition 2.2: 

Proposition 2.3. Assume that a solution to the discrete Schrodinger equation satisfies, Vj S 
7j, |uj(0)| + |uj(l)| < Ce““l , for some C > 0 and a > 0. Then, for t G (0,1) we have 
\uj(t)\ < Cie~°‘^ for all j G Z, with Ci not depending on t. 

Remark 2.1. When considering modified Bessel functions, the constant blows up at time t = 0 
and t = 1, and this is why first we have to prove the log-convexity in [toi^i] C (0,1) and then 
study what happens if to tends to 0 and ti tends to 1. In this case, the constant is independent 
of t and we can avoid this step in the justification. 


Proof of Proposition 2.3. Consider that Uj{t) = f{j,t). As in Proposition 2.1, thanks to the 
decay conditions we can extend f{x, 0) and f{x, 1) as entire functions and they are 27r—periodic. 
This means that at time t, f{z,t) is 27r—periodic and entire, and we recall that it is given by 


( 12 ) 


fiz.t) = f(x + iy,t) = 


g 22 i(cos a; cosh y —2 sin a: sinh y—1) 
g2z(i—1) (cos X cosh y —2 sin X sinh y—1) J 


Moreover, using Poisson’s summation formula 


|/(-z,0)|< ^ \uj{0)\e-^y Y 

j=: — CO j — — 00 

having the same estimate for \ f{z, 1)|. Now, if y > 0 we can write f{z,t) using the first and the 
second line in (fl^ in order to have that \f{z,t)\ < Cae^ Then, by Cauchy’s theorem, if 

j < 0, we have for all y > 0 that. 


Finally, we set y = —2aj, so \uj{t)\ < CaC af _ If J > 0 we can argue 

in the same fashion, changing the contour of integration. □ 


This proposition proves that under the hypotheses, the formal calculations are valid for Ga-e 
| uj(t)P, and by Fatou’s lemma we conclude the result. 


□ 
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In the three cases {[S, A]f, f) is written as the sum of two positive terms. We can use this 
fact and [HI (2.22)] in order to give a-priori estimates for solutions to (|T]). For example, when 
considering the weight , just by getting rid of the first sum in (fTTl) we have that 

(13) (cosh(4Ajfe) — dt < c(G(l) -I-G(O)), 

•’'O k=l 

and in particular this implies that in the interior we have more decay for the solution. On the 
other hand, using the formula 

d d 

k=l k=l 

cosh(4Ajfe)e^'^l^l \uj{t)\'^, 

k=l k=l 

we also have the bound 

(14) f t(l-t) ^ < c(G(l)-f G(0)). 

jGZ<i fc=l 

In the continuous case, in order to conclude Hardy’s uncertainty principle, this estimate for 
the gradient of the solution was crucial. However, here we can avoid the use of this estimate as 
we will see in Section 4, although it should be useful if one wants to relate this discrete result to 
the continuous one. Considering the modified Bessel functions, we can get similar estimates to 
those explained here for . Notice that in a £°° setting. Proposition 2.1 and Proposition 2.2 
imply extra decay in the interior of [0,1] for the solution as well. 


3. Log-convexity properties for solutions to perturbed discrete Schrodinger 

EQUATIONS 

When we introduce a potential in the discrete Schrodinger equation, we cannot use the method 
we use in the previous Section in order to justify the calculations. However, we can first prove a 
log-convexity property for solutions to @ where V = {Vj{t))j^id is a time-dependent bounded 
potential, and the solutions satisfy 

(15) + |wj(l)P) < +00- 

jGZ'^ 

Lemma 3.1. Assume that u is a solution to m where V is a time-dependent bounded potential. 
Then, for t G [0,1] and /3 G M we have 

jGZ'^ jGZ'i 

where C is independent of (3. 

Remark 3.1. ||l^||oo stands for supjgZ‘i,tG[o.i]{l^(OI}- 

Proof. We are going to assume, without loss of generality that /3 > 0. In order to give a rigorous 
proof of the result, we are going to truncate properly, following the procedure in m, the weight 
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SO that all the quantities that we are going to compute later on are valid and finite. To do 
this, we consider a function ip G such that 0 < </? < 1 and 


ip{x) 


1, a; < 1, 
0, x>2. 


Now, for iV S N we define (x) = tp (^) and 6^{s) = (3 (y))^dy, so 9^ G is 

non-decreasing and 




f3s, s < N, 
cP, s > 2N. 


Moreover, we have that |(0^)'(s)| < (3 and |(6*'^)"(s)| < Finally, we discretize 9^{s) by 
considering its evaluation at Z. In other words, = 0^(ji), for ji G Z. Notice that f Pji 

as N ^ oo. Now we take fj = e^^iUj{t) and compute the operators Sn and Aff, symmetric and 
skew-symmetric respectively such that dtfj = Sjvfj + Apffj + iVjfj. We have that 

{[Sn, AN]fp /,) = -23? ^ sinh(0j^+i - 29^ + 9l_^)f,+,J— 

jGZ'^ 

+ 2^2 (cosh(6»j^+i - sinh(6»j^+i - 9^) - cosh(6»^^ - sinh(6»j^ - |/jf, 


and we want to bound this quantity from below. To do that, we define ujv(a:) = 9^{x + l) — 9^{x) 
so that \92+i - 29^ + = |uAr(ji) - vnUi - 1)| < kw(C)l < l(6'^)"(6)l < ^ and 

sinh(6ij^_,_^ — 29^ + > — sinh • On the other hand, the factor that appears in the 

second sum is | cosh(i;Ar(ji)) sinh(wAr(ji)) — cosh(z;w(ji — 1)) sinh(uAr(ji —1)| < cosh(2^)|z;w(ji)— 
vnUi ~ 1)1 ^ since |^| < max{|p7v(ji — 1)|, kvOi)!} < /3- Combining these estimates 

we bound the commutator by 


{[SN,AN]fj, fj) > - sinh 


cp_ 

N 


+ ^ cosh(2,0) 


E I/. I 

jGZ'’* 


Now we use the perturbative version of Lemma 2.1 (see 0 ), taking into account that \dtfj — 
Sn/j - Ar,fj\ = \VM \< nioo|/,|. Thus, 

jGZ'^ jGZ'^ 

< gC(sinh(C/3/A)+/3cosh(2/3)/7V+||V|U) ^ (|u, (0) ^ + |u,(l)H- 

jGZ'i 


Using Fatou’s lemma, we conclude the result. □ 

Remark 3.2. Using the same method, it is straightforward to see now that, for A G R^*, 

(16) ^ E + l«.(l)l')’ 

jGZ'i jGZ'^ 


Once we have proved the lemma, it is quite easy to see that the log-convexity properties of 
Theorem 2.1, 2.2 and 2.3 are also satisfied when we add the potential V to the equation. Notice 
that this fact gives another proof of those theorems, just by setting the potential equal to 0, 
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although in this way we lose the a priori estimates piil) . (1141) . Let us see how we can do this in 
the case of Theorem 2.2. 

Theorem 3.1. Assume u = is a solution to the equation m where V is a time- 

dependent bounded potential. Then, for a > 0 and t € [0,1], 

n Yl {\uj {0)\'^+ \uj {1)1"^), 

fc=i jgZ'’* fc=i 

provided that the right-hand side is finite. 


Proof. When the right-hand side is finite, we have that (|mj( 0)P-I- |uj(l)p) is finite 

V Ai, A 2 G so applying (fTCll we have 

jGZ'^ jGZ'’* 

havingthat (Ai-|-A 2 )-j = X]fe(-^i,fc+^ 2 ,fe)-jfc- If we multiply this expression by 
and integrate it in (Ai, A 2 ) G the theorem holds using that 




□ 


In order to prove the same log-convexity properties for the other weights we can do the same, 
now using that 

Jr 

while, in order to get the inverse of the modified Bessel function we do not have an explicit 
formula, but it can be checked that multiplying the linear exponential by a similar function that 
we have used in the proof of Theorem 3.1 we get that the integral behaves asymptotically in 
the same way as the inverse of the modified Bessel function , whose asymptotic behavior is 
described in ([5]). 


4. Hardy’s uncertainty principle: Real variable approach 

In this Section we are going to give a discrete version of Hardy’s principle. As we have pointed 
out above, among the weights we have studied here we are going to follow the approach in 
considering the function . In order to do that we need the following Carleman inequality: 

Lemma 4.1. The inequality 

holds, for R large enough, when e > 0, > 0 and g G C'“(Z'^ x [0,1]). 
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Proof. Consider fj{t) = *)eiP (i+e)R^t(i gg^jj compute the oper¬ 

ators S and A symmetric and skew-symmetric respectively, such that 


Then, since 

ll^t/— *5/— ^/||^2([Q_i]_£2(2d)) > f {Stf + [S,A]f, f)dt, 

J 0 

what we need to prove is {Stf A [S,A]f,f) > ^-||/||f2(Z‘i)- After computing the commutator, 
we see that 

{Stf + [S,A]f, f) = sinh(2/r) ^ \ fj+et - fj-et P + 32^^ ^ fji + Rt{l - t) - |/jf 

j£Z‘‘ jGZ'^ ^ ^ ' 

-h2sinh(2^) ^ (cosh(4^(ji + Rt{l - t))) - l)\fj\‘^ - 32fA ^ (ji -H Rt{l - t))^|/jP 
jez’’' iGZ"^ 

+ 2fiR\l - 2tf Y. l/il" + - 2i) E cosh(2Ai(ji + 1/2 + Rt{l - t)))^{fj+eJj) 

jGZ'i jGZ'i 


sinh(2/r) EE \fj+ek fj-e 

jgZ<i k=2 


2 sinh(2/r) EE (cosh(4/rjfc) - l)\fj 

-,-gZ<i fc=2 


eR^ 

8 /r 


El/.- 


Thus, the commutator can be studied separately, according to the behavior in each variable 
jk, and the lemma holds if we prove that the following expression, written conveniently in one 
dimension, 

32^3 ^ ("j + Rt{l -t)- |/,f + 2^xR\l - 2t)2^ |/,f 

iGZ. ^ ^ ' j 

(17) -f 2 sinh(2/r) Y (cosh(4/r(j -h Rt{l - t))) - l) |/j ^ - 32/1^ Y 0 + “ ^)) Vi P 

3 3 

+ 8^^Ril - 2t) Y cosh(2Ai(j + 1/2 + Rt{l - t)))S(/,+i7‘) 
i 

is positive, since, going back to the expression for the commutator, this proves that the first 
three lines are positive, while the first two terms in the last line is clearly positive. By symmetry, 
the procedure will be the same if either t G [0, 5 ] or t G [ 5 , 1 ], so we assume t < 5 - What we 
are going to see is that when R is large, the leading terms in the expression for a general j make 
the sum positive. We have to distinguish several regions and possibilities in order to give the 
whole proof of the lemma. First, let us assume that t > ^ — 3^^23727 M > Ao, where /io is the 

positive value such that cosh^(/xo) = 2 sinh( 2 ^o )7 having that /xq — 0.255413_In this case we 

use 2Q{fj+ifj) > —l/j+ip — l/jp to bound the expression by 

82^^^Y (j+Rti^ -t)- I/.I' + 2mA 2(1 - 2t)"E l/.l' 

-h 4sinh(2^)^ (cosh2(2/x(j -f Rt{l - t))) - l)|/jp - 32^^^ (j + Rt{l - t)f\fj\^ 

3 3 

— 8 /iCOsh(/x)i?(l — 2t) cosh(2^(j + Rt(\ — t)))|/jP, 

3 
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SO in this way we have an expression of the form -F(/r, R, and we would like to have that 

F is positive. If j is given such that cosh(2^(j + Rt{l — t))) > then, the leading 

terms of F, when R is large, are 

4sinh(2^) cosh^(2^(j + Rt{l — t))) — 8/iCOsh(/i)i?(l — 2t) cosh(2/r(j + Rt{l — t))) 

> c„sh(2Mi + «(1 - tm ( 5 ^ - = 0. 

On the other hand, if j is given such that cosh(2/x(j + Rt{l — t))) < then, since 

this implies that j + Rt{l — t) is much less than we see that the sign of F is determined by 
^2 

-8/iCOsh(/r)i?(l — 2t) cosh(2^(j + Rt(l — t))) > 0. 

8 ^ 

Hence, we have to consider now j such that cosh(2/i,(j + Rt{l — t))) = aR, with a between 
the two values of the previous conditions. But in this case the leading terms in the expression 
are written as 

+ 4sinh(2Ai)a2 - , 

and this is a polynomial of degree two in a, whose discriminant is negative since /i > /loi so 
there is no change of sign in the polynomial and it is positive. This completes the proof when 
16 ^ ^ > ^ 0 - 

When t < ^ — J M > Mo, we use the same reasoning to show that if j is given such that 

(18) 

cosh(2^(j + i?t(l-<))) < or cosh(2/r(j + i?t(l-t))) > 

then the expression F is positive for R large enough. On the other hand, when j is not in that 
region, we cannot use the same argument as before, because now we can have negative terms. 
However, we can define jo and Jq such that 

• jo+Rt{^—t) is positive and for the first time cosh(2/r(jo+i?t(l—<))) > g4^2^cosh 

• jg+i?t(l—t) is negative and for the first time cosh(2^(jg+i?t(l—<))) > 2t). 

We are going to see here that for j = jo we can absorb the negative terms in CZD- Furthermore, 
we can use the same argument to prove the same for j = jo + 1, jo + 2 ,... until we go into the 
region defined by (fTHll . Actually, it is easy to check that jo + 3 is in the region defined in (fT^ . 
Moreover, the proof for j = jo, jo — 1,.. ■ is the same. 

Since jo — 1 is in m we have that the part of CZD involving fj^ that we need to study is 
(19) 

|^32/r3 (^jo + Rtil -t)- + 2f,R^l - 2tf - Z2y.\jo + Rt{l - t))'^ \fj‘^ 

+ (4sinh(2/r)( cosh^(2^(jo + Rt{l — t))) — l) — 4/ii?(l — 2t) cosh(2/r(jo — 1/2 + Rt{l — t)))) |/joP 
+ 8^iR{l - 2t) cosh(2/r(jo + 1/2 + Rt{l - t)))^{fjo+ifjo) 

+ 2sinh(2/x)(cosh^(2^(jo + 1 + Rt{l - t))) - l)|/jo+ip. 
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where we add a part of cm involving fjg+i in order to gain an extra positive term. Now first 
we see that if ^{fjo+ifjo) — 0) then (fTUl) is positive. Indeed, in this case we get rid of the last 
two lines in m, and considering that cosh(2/r(jo + Rt{l — t))) = aR{l — 2t) we can write the 
leading term in m as 

2fiR^{l — 2t)^ + 4sinh(2^)a^i?^(l — 2t)^ — 4/io; ^ —, 

cosh fi + sinh fi 

since cosh(a; + a) = cosh x cosh a + Vcosh^ x — 1 sinh a ~ cosha::(cosha + sinha) => cosh(a; —a) ~ 
cosh a^sinh a ' Computing the discriminant of this polynomial we see that the expression is 

positive. 

On the other hand, if ^{fjg+ifjg) < 0, there is zq G C such that 5(zo) = 7 o > 0 and 
fjo+i = -zofjo ^(/jo+i/jo) = -7o|/,oP: l/jo+iP > 7ol/ioP- By using this we can bound the 
leading term in (dH) by 

— + 2fiR^{l — 2t)‘^ + 4sinh(2^)a^i?^(l — 2t)^ — ^ — 

8 /r cosh /r + sinh fi 

— 8^R^{1 — 2t)^a(cosh /i + sinh/r) 7 o + 2 sinh(2^)a^i?^(l — 2t)^(cosh(2/i) + sinh( 2 ^))^ 7 Q. 


Since 1 > 1 — 2t we can write this as R^{1 — 2tY times a polynomial expression in a and 7 
of degree 2 in each parameter. We can write it as a polynomial of degree two in a and compute 
minus the discriminant, getting 


( 20 ) 


4 



^ (4sinh(2^) + 2 sinh( 2 ^)(cosh( 2 /r) + sinh( 2 ^))^ 7 g) 


^ 8 /i 7 o(cosh p, + sinh /r) 


cosh ^ + sinh fi 


2 


We want the latter to be positive, and this is another polynomial, now in 70 , of degree two, 
so we compute again the discriminant of the polynomial and we get 

4096/i"^ — 4 f32^sinh(2^) — --^^^ ^ sinh( 2 /x) \ /-p sinh( 2 /r))^ 

\ (cosh p, + sinh/j,)^ /j, J \ 

+ (cosh(2/x) + sinh(2^))^ — 64pt^(cosh/r + sinh/i)^^ , 

M / 

which is easy to see that is negative. Hence, we have that (EHl) is positive and therefore the part 
of (1171) we are studying is positive. As we have pointed out above, we can repeat this argument 
for the rest of coefficients that we have to check to see that (El) is positive in this case. 

It only remains the study of the case /r < /tq- In order to proof Theorem 4.1 below, we do not 
need this case, since we are going to assume fi > 1/2 > Nevertheless, it is easy to see that 
we can follow the following sketch to see that in this case (1171) is also positive. 

• Assume j such that cosh(2/r(j + Rt{l — t))) < After identifying the leading terms 
in (ED for that j, we get R^ times a polynomial of degree two in the variable (1 — 2t), 
whose discriminant is negative when /i^ cosh^ /i < 625. Since fi < fiQ we are in this case 
and the negativity of the discriminant implies that the leading term is positive. 

• If cosh(2/x(j + Rt{l — t))) > R then we can prove the positivity of the leading term using 
4sinh(2/r) > 8fi cosh fj,. 
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• Finally, if cosh(2^(j + Rt{l — t))) = aR with a G [^, l], the leading term leads again 
to a polynomial expression of degree two, which is negative when 

64/1^ cosh^ /I — 16 sinh(2/r) ( — + 2/r ) < 0, 

and this happens when fj, < ^q. Again, the latter implies that the leading term is positive. 

□ 


After proving this Carleman inequality, we can give a version of Hardy’s uncertainty principle 
in the discrete setting using the weight as the discrete version of Gaussian decay, following 
the procedure used in [S]. 

Theorem 4.1. Assume u = is a solution to the equation where V is a time 

dependent bounded potential, which satisfies the decay condition (0, for some A > 1/2. Then 
u = 0. 


Proof. As we have seen in Section 3, the decay condition (fTIl implies that there is N\y depending 
only on A and ||H||oo such that 


sup ^ < Nxy 

‘6 [ 0 . 1 ] 


For given i? > 0, we take pL and e such that 


(l + e)3/2 
2 ( 1 -e)3 


< p, < 


A 

1 + e’ 


and we consider, for i? < M G N, smooth functions 0 < 9^ (x), r]ii{t) < 1 such that 


e^{x) 


1, |x|<M, 

0 , |a:| > 2M. 


, Vnit) 


1 t e [A 1 _ A] 

0, t G [o, U [l — l] . 


Then we define 6^ = 9^(j) for j G and gj{t) = 9^r]ji{t)uj{t), so gj{t) is compactly 
supported in Z'^ x (0,1) and satisfies 

d 

dtgj - i{Adgj + V^gj) = PROfufit) - - {9f^ - 
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We apply Lemma 4.1 to show that there is a constant W such that 

( 21 ) 

Jo ^ ~ ^ Jo , 

Jo ^ 

«1 rf 

+ ^ / ^^g2pb-efc+Rt(l-t)ei|"-(l+e)fl"t(l-«)/8/^^2^|0M _ HujP 

'■^0 j fe=l 

-'o „ 


i fe=i 


Now we have to understand each term of the right-hand side of m separately. The first 
term can be hidden in the left-hand side by taking R large enough. For the second term, it is 
supported where —2M < j < 2M and t G [^, U [l — 1 — j so, since 2ab < ea? + fo^/e, 

mIj + -^^(1 ~ ^)eiP ^ m( 1 + e)|jP ^ g) “ 7’ 


while, for the other two regions, the support in t and j changes, but the same reasoning leads to 

^Aj ~ + Rt{l — t)eip < A|jp H— {Rt{l — t) + 1)^ < A|jp -I-, 

e e 

^k\j + Cfe -f Rt{l — t)ei p < A|j p H — {Rt{l — t) + 1)^ < A|jp -I-, 

e e 

with c independent of R, A, and e. Thus, using the natural bounds for the functions 9^ and ry^j 
we get 


/ E' 




( 22 ) 


< W(i?^e2^/'=-b sup V 

t 6 [ 0 . 1 ]“ 


Furthermore, the term multiplied by M~‘^ tends to zero when M tends to infinity. On the 
other hand, we can bound the left-hand side just integrating in t G [(1 — e)/2, (1 + e)/2] and 
considering that |jj < e(l — e)^i?/4. Under these circumstances, gj = Uj and moreover 

g\j + Rt{l - t)eip - (1 -be)i?^t(l - t)/16g > - ef - (1 + e)^) > 0, 

so there is C{X, e) > 0 and N\^v,e only depending on A, ||U||oo and e such that 


„C(A,£)R= 


ll“llL=((i^,i^),^2(b-|<£(l-£)2R/4)) - 


Combining the fact that N ^||m(0)||£2 < ||u(f )||^2 < iV||u( 0)|72 for 0 < t < 1 and N = 
with 

l|u(t)|U2(z<i) < ||u(t)|U2(|j|<e(i_g)2fl./4) +e~^'' ^ ^^^Nxy, 

we conclude that 


|w( 0 )ll£ 2 (z<i) < Nxy^^{e 


-C{\,e)R^ I -Ae^(l-e)*R' 


/16)^ 
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Finally, we let R tend to infinity to have that u = 0. □ 

According to (O and [T^, this is not the best result we should expect. In the one dimensional 
setting, the sharp results assume decay conditions of the type with a<l,tiG{0 ,l}, 

so it is reasonable to think that in the multidimensional setting the sharp condition is similar, 
now changing the modified Bessel function Ij{a) with the product of modified Bessel functions 

11^=1 ^ 3 k (“)■ 
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